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We study theoretically the electronic entropy and specific heat in diffusive superconductor-normal 
metal-superconductor (SNS) Josephson junctions. In particular, we consider the influence of non- 
idealities occurring in an actual experiment, such as the presence of barriers at the NS interfaces, 
the spin-flip and inelastic scattering in the N region and quasiparticle subgap states in the supercon- 
ductors. We find that spin-flip and inelastic scattering do not have, for typical parameters values, a 
large effect. On the contrary, the presence of barriers suppresses the superconducting correlations in 
the N region, with the consequence that the entropy and the specific heat get reduced eventually to 
those in the absence of superconductivity for opaque interfaces. Finally we suggest an experiment 
and check that it is possible, under realistic conditions, to measure the dependence of electronic 
specific heat on the phase difference between the superconductors. 
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I. INTRODUCTION 

When a superconductor (S) is put in contact to a nor- 
mal metal (N), Cooper pairs can leak across the inter- 
face. As a result there exists a non-vanishing pair am- 
plitude in the normal metal. This is the so called prox- 
imity effect at the heart of the properties of SN hybrid 
systemsii Already studied during the early days of super- 
conductivity, the interest in proximity devices has been 
renewed when it became possible to study hybrid systems 
with dimensions comparable or smaller than the coher- 
ence length in the normal metal. Over the last years 
an intense theoretical and experimental activity has un- 
veiled numerous interesting aspects of proximity-based 
devices (see RefsiUliil^- So far most of the attention 
has been devoted to their transport properties, although 
their thermal properties are of great importance espe- 
cially in relation to the use of hybrid devices for electron 
coolingi^ Thermodynamic properties have been consid- 
ered in Ref]^ for a SNS structure showing that proximity 
effect may lead to a substantial deviation of the specific 
heat from that in the normal state and that it can be 
largely tuned in magnitude by changing the phase dif- 
ference between the superconductors. To confirm these 
predictions a measurement setup, which makes use of a 
highly-sensitive recently developed method^^iii has been 
also suggested. 

In the present paper we further elaborate on this prob- 
lem. Namely, we extend the results reported in Refl^ in 
order to investigate the effect of non-idealities occurring 
in an actual experiment, such as inelastic and spin- flip 
scattering, and the role of the barriers and the SN inter- 
faces. Moreover, as most probably the experiment will 
be performed embedding the SNS junction within super- 
conducting loops, the signal will also collect the contri- 
bution due to the superconducting region, as well as the 
unavoidable contribution due to phonons. We include 
such additional contributions in order to find the condi- 



tions under which they do not hinder the signature of 
proximity effect. 

The paper is organized as follows. In Sec. |TT] wc de- 
scribe the general formalism based on the quasiclassical 
theory for diffusive SNS junctions. In Sec.|IIT]we present 
our results. We show the entropy and specific heat of 
the proximized N region in a SNS junction calculated as 
functions of temperature and superconducting phase dif- 
ference, and compare with the ones relative to a metal 
in the normal and superconducting state. Furthermore 
we explore the dependence of the entropy and specific 
heat on the strength of spin-flip and inelastic scattering, 
interface resistance and presence of subgap quasiparti- 
cle states. Moreover, we check that the proposed ex- 
perimental setup allows to measure the phase-dependent 
contribution to the specific heat of the SNS junction. We 
summarize our main conclusions in Sec. IIVI 



II. THEORETICAL FRAMEWORK 

In this section we will introduce the theoretical for- 
malism needed to calculate the entropy and the specific 
heat of the system depicted in Fig. [1] which consists 
of a superconducting loop interrupted by a SNS junc- 
tion. The N metal length is equal to L, and the phase 
difference cf) across such junction is controlled through a 
perpendicular magnetic field H . 

The formalism of quasiclassical Green's functions has 
been proved to be a powerful tool to describe mesoscopic 
superconductivity In the dirty limit the mean free path 
is much smaller than coherence length and the Eilen- 
berger equations can be simplified to the Usadel equa- 
tions^ which we use in this paper. Although much sim- 
pler than the original microscopic version, the equations 
for quasiclassical Green's functions are still hard to solve 
analytically and numerical methods are usually unavoid- 
able. The following two different parametrizations for 
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FIG. 1; (Color online) Schematic of a SNS junction embedded 
in a superconducting loop. The normal (N) metal region, of 
length L and width W, is placed between two barriers and 
subject to inelastic and spin-flip scattering whose rates are, 
respectively. Fin and F^f. Fs is the inelastic scattering rate in 
the superconductor (S). The perpendicular magnetic field H 
and the geometry allow a superconducting phase difference (p 
to appear between the S electrodes across the N region. $ 
is the total flux threading the loop area and x = denotes 
the middle of the junction. The length values indicated are 
typical for a system made of a Cu as N and Al as S materials. 



where L is length of the normal region, A is the 
temperature-dependent superconducting order parame- 
ter, and Fs takes phenomenologically into account a pos- 
sible smearing of the density of states (DOS) due to 
inelastic scattering in the superconductor The latter 
might be due either to impurities or to inverse proxim- 
ity effect induced by the nearby-contacted metal.— We 
recall that (j) denotes the phase difference between the 
S electrodes. For the purpose of the present paper, the 
self-consistent solution of the order parameter is not en- 
forced so that A is assumed constant in S and zero in 
the N region. Furthermore, we assume the usual BCS 
temperature dependence for A with critical temperature 
Tc = Ao/(1.764fcB), where Ao is the zero-temperature 
order parameter, and fee is the Boltzmann constant. If 
_Rn is the resistance of the N region coupled via a barrier 
with resistance i?B to a superconductor, the boundary 
condition for 9, derived by Kupriyanov and Lukichev,— 
becomes 



rLd.JU=±L/2 = ±smh[e{±L/2) - 9s] 



(3) 



where r = R^/R^, and 6s is defined in Eq. The 
local DOS in the normal region can be obtained by 

N'ix, E, T, (f) ^ AApRc {cosh [e{x, E, T, 0)]} , (4) 

where TVf is the DOS per spin at the Fermi level in the 
absence of superconductivity. 



the quasiclassical equations arc often usedi^ The first 
one, and most widely employed, is the 9 paramctriza- 
tion which, however, turns out to yield unstable differ- 
ential equations in certain cases (for example when the 
phase difference is close to tt). The second one is the 
Riccati parametrization which is particularly suited for 
time-dependent problems, and yields stable differential 
equations. In the next subsections we will briefly intro- 
duce these parametrizations, and define the electronic 
entropy and the specific heat of the system. 



A. 9 parametrization 

In the 9 parametrization^ the Usadel equations can be 
written as 

hDdl9 = •-2i{E + iFin) sinh(6l) + ^ {d^xf sinh(26l) 
sm\i{2e)dJd^X + &\nh^{9)dlx = 0, (1) 

where 9 and x ^-re complex scalar functions of energy 
E and position x, D is the diffusion coefficient and Fjn 
describes a small inelastic scattering rate in the N region. 
The energy E is measured from the chemical potential of 
the superconductors. The boundary conditions at the NS 
contacts for perfectly transmitting interfaces are^^ 



e{±L/2) = 6*3 = arctanh[A/(£; + iFs) 
x(±L/2) = ±0/2, 



(2) 



B. Riccati parametrization 

The Riccati paramctrizatio n"'^^'^^ has been employed 
by Hammer et ali^ in order to calculate the supercurrent 
as well as the DOS for SNS junctions in the presence of 
opaque NS interfaces and spin-flip scattering. Here we 
just present the main formulas for completeness. The 
Usadel equations read 



+2i(£: + iFi„)ai,2 = 0, 



where 



± 



I — aia2 
1 + aia2 



2a2 



1 + Q!lQ!2 



(5) 



(6) 



and ai and a2 are complex scalar functions of energy 
E and position x. Here we have also introduced an ef- 
fective spin-flip scattering rate Fgf which accounts for 
pair-breaking mechanisms, such as due to magnetic im- 
purities or to the presence of a perpendicular magnetic 
field H. For a thin normal region with small width 
W (smaller than or comparable to the coherence length 
in the N metal ^ = y^hD/ A), such rate is given by 
Fst = De'^H'^W^/{6h),.i^ The boundary conditions for 
Eqs. jl]), which can include the transmission coefficient of 
the interface and the ratio r, are given in Refill. Finally, 
the local DOS is given by 



JV{x, E, T, (j)) = TVpRe [^i {x, E, T, i 



(7) 
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C. Electronic entropy and specific heat 

The electronic specific heat (C) of the proximized N 
region of the SNS junction in Fig. [1] can be calculated 
through the expression 



C(T,, 



T 



dSjT, , 



(8) 



where the electronic entropy S is defined as 

S{T,(f>) = f- / dxdEAf{x,E,T,(l)) (9) 

^ J-L/2J0 

X {fiE) Hf{E)] + [1 - f{E)] ln[l - f{E)]} , 

and f{E) = {1 + exp[£;/(fcBT)]}"^ is the Fermi-Dirac 
quasiparticle distribution function at temperature T. We 
introduce the position-dependent specific heat (c) and 
entropy (s) according to the expressions 



5(r,, 



L J-L/2 
1 



L/2 



dx c{x, T, ( 
dx s[x, T, ( 



(10) 



In order to allow a comparison of our results with the 
experiments it is important to calculate the total spe- 
cific heat of the loop. To this end one should also add 
the electronic contribution of the superconducting part 
of the ring (which can be rather small due to the ex- 
ponential suppression with temperature), as well as the 
contributions from the lattice phonons of the N and S 
parts. We will discuss this issue at the end of the paper. 



III. RESULTS 

In this section we report our numerical results for en- 
tropy and specific heat determined by using Eqs. ([5]) and 
^ and the expression for the DOS. The latter is calcu- 
lated by numerically solving the Usadel equations, using 
the relaxation method^ in one of the above parametriza- 
tions. 

In the following analysis we shall focus on long SNS 
junctions only (Ag 3> Ej-h): a-s this is the relevant regime 
for metallic diffusive Josephson junctions. The main gen- 
eral feature is that the DOS is an even function of energy, 
exhibiting a minigap Eg^ and it is position-dependent. 
Moreover, Eg depends on Thouless energy, phase dif- 
ference, inelastic and spin-flip scattering rate, and ra- 
tio rJ^ In particular, Eg is maximum at (/> = 0, where 
Eg « 3.12i?Th in the limit of long SNS junctions, whereas 
it vanishes at (f> = tt ^ As we shall see, entropy, and there- 
fore specific heat, are heavily affected by the magnitude 
of the minigap, which is the most evident manifestation 
of the proximity effect. 

The electronic entropy of the proximized N region of 
the SNS junction versus temperature is compared with 
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FIG. 2: (Color online) (a) Entropy <S vs temperature T com- 
pared to that of a normal metal (N) and of a superconductor 
(S) calculated for Ao = 10i?Th and (/> = 0. (b) Local entropy 
s vs T at different positions x along the N region calculated 
for Ao = lO^Th and = 0. We have set r = 0, Fin = 0, 
Fsf = and Fs = 0, while S and s are normalized with re- 
spect to the N metal value. Tc is the critical temperature of 
the superconductor. 



that of a normal metal (N) and of a superconductor (S) 
in Fig.[2{a). Here we set Aq = IOStii and = and nor- 
malize S to the N case value, namely iSn = 7T, where 7 = 
27r^A/'Ffc|/3. For the case of a superconductor, the figure 
shows the standard BCS result^ where the electronic 
entropy is exponentially damped for T <^ Tc, stemming 
from the presence of the gap A in the DOS. Figure UK a) 
shows that in the proximized region S strongly departs 
from the normal-metal behaviour especially at low tem- 
perature. In particular, for T <C Tc, 5 exhibits a strong 
suppression that resembles that one of a superconductor. 
It turns out that the temperature determining the onset 
of such suppression is set by the proximity-induced mini- 
gap Eg. More precisely, for k^T <^ E^i-^ and in the long 
junction regime (Aq 3> i^Th), for = the electronic 
entropy can be described rather well by the following ex- 
pression 



5(T, = 0)= 5Nexp(-a.BTh/fcBT), 



(11) 



where for the prefactor we get a « 2.1. Accordingly, 
in the same limits we get the following approximated 
expression for the electronic specific heat: 

C(r, = 0) ~ CN(aSTh/fcBT)exp(-aSTh/fcBT), (12) 

where Cn = jT is the specific heat in the normal state. 

In Fig. the position-dependent electronic entropy 
s{x,T, 0) is plotted as a function of temperature for dif- 
ferent values of x. First of all we note that at the NS 
interface {x = L/2) s coincides with that of the super- 
conducting lead. Secondly, by moving toward the centre 
of the N region, s tends to approach the entropy of a N 
metal, while the temperature which determines the onset 



4 




TIT^ x/L 



FIG. 3: (Color online) (a) Entropy <S vs T at = calculated 
for several ratios AoZ-Exh- (b) s vs x at = and T = 
0.5-ETh/fcB calculated for several ratios AqZ-Etii- We have 
set r = 0, Fin = 0, V^i = and Fs = 0, while S and s 
are normalized as in Fig. [2] N and S curves are reported for 
comparison. 




FIG. 4: (Color online) (a) Entropy <S vs T calculated for 
Ao = lO-ETh and different values of (j> . (b) Entropy iS vs 
calculated for Ao = IO-Etii and several values of temperature. 
We have set r = 0, Fin = 0, Fgf — and Fs = 0, while S is 
normalized as in Fig. [2] N and S curves are reported for 
comparison. 



of strong suppression gets reduced. This latter behaviour 
stems from the weakening of the proximity effect (i.e., re- 
duction of superconducting correlations) departing from 
the NS interfaces. 

Figure [3]^a) shows the value of S versus T calculated 
for several lengths of the N region expressed by the ra- 
tio Ao/-Etii at = 0. As expected, for long junctions 
(Ao = lOOE'Th) the curve approaches the N metal en- 
tropy, which follows from weakening of the proximity ef- 
fect in longer SNS junctions. Similarly, also the tem- 
perature which sets the onset of strong suppression gets 
reduced for longer junctions, as the minigap in the DOS 
gets suppressed by increasing L. In the opposite limit 
(Aq = 5£^Th), the entropy approaches that in the S state 
so that shorter SNS junctions are more suitable in order 
to maximize the effect of proximization on the entropy. 

The behaviour of the position-dependent electronic en- 
tropy s for fixed length L allows to investigate the role 
of different superconducting gaps, and it is shown in 
Fig. Eljb). Here we set T = 0.5£;Th/fcB and </) = 0. The 
entropy strongly decreases by moving away from the mid- 
dle of the junction (x = 0), approaching the value in the 
superconductor at the NS interface [x = L/2). We em- 
phasize that the increase of the value of Aq leads to an 
overall strong suppression of the entropy. This can be 
understood by noting that, for a given L, smaller-gap su- 
perconductors reduce the proximity effect in the N region 
producing an enhancement of the electronic entropy. 

Figure m^a) displays the electronic entropy iS as a func- 
tion of T for Aq = 10i?Tii and several different values of 
(j). The curves for = tt were obtained through the 
Riccati parametrization which allows stable solutions of 
the Usadel equations. The data show a strong phase- 
dependence of S stemming from the coherent nature of 



proximity effect. In particular, by increasing (j) from zero 
one finds an entropy enhancement in the whole temper- 
ature range mainly due to the suppression of the mini- 
gap Eg. Remarkably, at = tt the minigap vanishes 
and, for temperatures below T ~ 0.15Tc, the entropy 
shows a linear behaviour, namely 1S/7 = 0.65T, which 
corresponds to a 35% reduction with respect to the en- 
tropy of a N metal S^. The full phase dependence is 
shown in Fig. Sfb), where S versus phase is plotted for 
Ao = 10i?Tii at several different temperatures. In par- 
ticular, for each temperature the electronic entropy in- 
creases by increasing the phase difference between the 
superconductors, and it turns out to be maximized at 
(j) = TT, i.e., when the minigap is zero. We stress that the 
overall entropy variation with (f>, especially at the lowest 
temperatures, can be very large. This directly reflects on 
the electronic specific heat phase dependence. 

Figure[n{a) shows the electronic specific heat C vs tem- 
perature for a junction with Ao = 10i?Th at different val- 
ues of the phase difference between the superconductors. 
In this calculations we extended the results of RefH up 
to = TT. A strong phase-dependence is found at all 
temperatures: by increasing (f> from zero to tt, C tends to 
increase for small temperatures and to decrease for tem- 
peratures close to Tc- A linear temperature dependence 
is found for (/) = TT below T = 0.15Tc. 

In the next subsections we will discuss the influence of 
inelastic and spin-flip scattering, the presence of a barrier 
at the NS contacts, and the presence of subgap states 
within the gap induced by a finite value of Fs. For the 
sake of definiteness, in the following we will focus on the 
case = and Aq = 10i?Th, which shows the largest 
deviation from the N metal case. 



5 




FIG. 5: (Color online) (a) Electronic specific heat C vs T cal- 
culated for Ao = 10-ETh and different values of 0. (b) C vs 
calculated for Ao = lO-Exh and several values of temperature. 
We have set r = 0, Tin = 0, Tsf = and Ts = 0, while C is 
normalized with respect to the N metal value. N and S curves 
are reported for comparison. 



A. Inelastic scattering 

Inelastic scattering of quasiparticles in the N region is 
introduced through the parameter Fin appearing in the 
Usadel equation. Here we assume that, at low tempera- 
tures, the inelastic scattering is independent on tempera- 
ture. A non-zero value of Fin affects the thermodynamic 
quantities by inducing a finite DOS below the mingap. 
Figuresinja) and[6l^b) show the entropy and specific heat, 
respectively, as a function of temperature for several dif- 
ferent values of Fin. The curves corresponding to the 
normal and superconducting states are displayed for com- 
parison in the same plots. As expected, for increasing 
inelastic scattering the entropy increases at all tempera- 
tures, while the specific heat increases for T < 0.2Tc but 
decreases for larger temperatures. Note, however, that at 
low temperatures Fin is typicallyi^ of the order of 10~'^A, 
which corresponds to Fin ^ 0.01i?Th when Aq = 10i?Th- 
In turns out that such small value has negligible effect on 
the entropy and specific heat of the system. 



B. Spin-flip scattering 

We consider here the effect of pair-breaking mecha- 
nisms by setting a finite spin-flip rate Fgf in the normal 
metal, while assuming that Fgf = in the superconduc- 
tors. According to Reffisl. in the presence of spin-flip 
scattering the DOS is suppressed and the minigap gets 
reduced, disappearing eventually. One should, therefore, 
expect a change in the entropy and specific heat. Fig- 
ures [Tlja) and (b) show the plots of S and C, respectively, 
for various values of Fgf . Indeed, by increasing Fgf , the 
entropy increases for all temperatures while the specific 
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FIG. 6: (Color online) Influence of inelastic scattering on (a) 
entropy S and (b) specific heat C for Ao — IOBtii, = and 
several values of Fin. We have set r — 0, Fgf — and Fs = 0, 
while S is normalized as in Fig. [2] and C as in Fig. [S] N and 
S curves are reported for comparison. 



heat increases only in a smaller temperature range. More 
precisely, at T « 0.2irc the spin-fiip scattering does not 
modify the specific heat as all plots cross each other in 
this point [see Fig. [7{b)]. We remark that, as stated 
in Sec. Ill Bl the spin-flip rate is related to the intensity 
of the magnetic field applied perpendicularly to the loop 
and needed to establish the phase bias. Using the nu- 
merical values given above we find that Fgf equals i?Th 
only for a large magnetic field. Namely, Fgf = Etu for 
= 281G = 2.6Hc, where He = 105G is the aluminium 
critical magnetic field. For small magnetic field no rele- 
vant changes must be expected. 
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FIG. 7: (Color online) Influence of spin-flip scattering on (a) 
entropy S and (b) specific heat C for Ao — 10i?Th, = and 
several values of Fsf. We have set r = 0, Fin = and Fg = 0, 
while S is normalized as in Fig. [2] and C as in Fig. [5] N and 
S curves are reported for comparison. 
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FIG. 8: (Color online) Influence of non-ideal interfaces on (a) 
entropy S and (b) specific heat C for Ao = lO-Exh, = and 
several values of r. We have set Fin = 0, Fgf = and Fs = 0, 
while S is normalized as in Fig. [2] and C as in Fig. [S] N and 
S curves are reported for comparison. 




FIG. 9: (Color online) Influence of quasiparticle subgap states 
on (a) entropy S and (b) specific heat C for Ao = lO-Exh, 
(j> = Q and several values of Fs = 0. We have set Fin = 0, 
Fsf = and r = 0, while <S is normalized as in Fig. [2] and C 
as in Fig. [5] N and S curves are reported for comparison. 



Non-ideal NS interfaces 



E. Contribution of the superconducting loop 



The presence of a barrier at the NS boundaries is taken 
into account through the ratio of resistances r and the 
transmission probability of the barrier r. Their effect, 
on both DOS and supercurrent, in a SNS structure has 
been considered in Reffisl. It was found that, on the 
one hand, by increasing r from zero not only the DOS 
but also the minigap gets suppressed. On the other, the 
minigap only slightly depends on the value of r. This 
motivates us to explore the effect of a finite r on entropy 
and specific heat, while keeping r = 1. In Figs.[5{a) and 
m^b) we report the plots of entropy and specific heat, 
respectively, for various values of the ratio r. It can be 
noted that both S and C are strongly affected by a finite 
value of r. Indeed, already for r = 3 the curves are 
quite close to the one relative to N apart from very small 
temperatures. We conclude that very small values of the 
ratio r are necessary to observe a clear phase-dependent 
specific heat. 



D. Subgap states 

The effect of the presence of quasiparticles states 
within the gap can be studied by setting Fs 7^ in 
Eq. ([2]). Figures HJa) and [HJb) show the plots of en- 
tropy and specific heat, respectively, function of 
temperature for various values of Fg. The effect of a 
finite Fs is very weak up to Fg « O.li^Th- Since the typ- 
ical experimental values are between Fg = 10~'^Ao and 
Fg = 10~^Aor^ these results demonstrate that virtually 
no influence can be expected from quasiparticles states 
within the gap. 



Electronic specific heat can be measured by ac 
calorimetry as discussed in Refl24l Recent developments 
in highly-sensitive calorimetry applied to mesoscopic sys- 
tems allow to measure small signals J^iii In this section 
we estimate the specific heat which would be measured in 
an actual experimental setup, including the phonons con- 
tribution, aimed at observing the tuning of the specific 
heat by changing the phase difference between the S elec- 
trodes. To determine the total specific heat we suppose 
that the system consists of a SQUID loop of square shape 
(i.e., like the one depicted in Fig.[T]) with side length equal 
to 1/im, N region length L ~ 0.6/xm, width W ~ 0.1/im 
and constant thickness. We assume that the normal ma- 
terial is Cu while the superconductor is Al, and that the 
superconducting coherence length ^ ~ 190 nm. The total 
specific heat C"°* can be calculated through the expres- 



Al I ,.Alrp3\TAAl 



where the second terms in the round brackets account for 
the low-temperature phononic contribution. Here 



,.Cu(Al) 



27r^ 



5 (fiwCu(Al))3^ 



(14) 



^Cu(Ai) lacing the sound velocity in Cu (Al)i^ F*^" and 
V^^ are the volumes relative to the N region (Cu) and 
the S region (Al), respectively. is the total vol- 

ume of the structure. In Figs. [TO] we collect the re- 
sults for the total specific heat normalized to its nor- 
mal state value calculated at Al critical temperature Tq: 

(^tot^tot ^ (^Cuj.^ _^ ^Cuj.3)yCu ^ (^Aly^ _^ ^Alj.3)yAl 

for Ao = 10i?Tii- As shown in Fig. [TOKa). where plots of 
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FIG. 10: (Color online) Total specific heat of the loop (7*°' , 
calculated including the contribution from the S regions and 
the phonons, as a function of (a) temperature and (b) phase- 
difference. C'°* is normalized with respect of the normal-state 
total specific heat evaluated at the critical temperature of Al. 
Plots are relative to the following parameters: Ao = IOBtii 
and r = 0. 

C*°* as a function of temperature are reported for various 
values of the overall specific heat maintains a differ- 
ent functional dependence on T for the different values 
of phase-difference. In Fig. [TUlJb) C*°* is plotted as a 
function of for various temperatures. The comparison 
between Fig. [IW b) and Fig. [5]Jb) makes clear that the to- 
tal specific heat is still strongly phase-dependent, despite 
the inclusion of the contributions from the S regions and 
the phonons. We conclude that phase-dependent specific 
heat can be measured, and that can be maximized in the 
case of N regions of intermediate length (Aq = lOSxh)- 
We emphasize that while in the case of very long junc- 
tions (Aq ^ -ETh) the specific heat is close to the normal 
state value and weakly phase-dependent, for short junc- 
tions (Ao ^ i?Th) the specific heat of the S regions dom- 
inates therefore masking the dependence on the phase. 



perature and superconducting phase difference aiming 
at investigating the role of non-idealities occurring in an 
actual experiment. Indeed, the presence of barriers at the 
NS contacts, as well as finite inelastic and spin-flip scat- 
tering produce a weakening of the proximity effect which 
is expected to affect both entropy and specific heat and 
to suppress their phase-dependence. 

Here we have focused on the intermediate junction 
length, with Aq ~ lOi^Tiu which turns out to be a fa- 
vorable condition as far as phase-dependence visibility is 
concerned. As a general feature, we have shown that, 
at low temperatures, the entropy presents an exponen- 
tial suppression which resembles the one of a supercon- 
ductor and where the role of the superconducting gap 
is played by the minigap. The most important results 
arc the following. First, we have found that inelastic 
scattering in the N region negligibly affect the specific 
heat, at least at the typical experimental values relative 
to low temperatures. We arrive at the same conclusion 
also for subgap quasiparticle states in the superconduc- 
tors. Similarly, the presence of spin- flip scattering, re- 
lated to the applied perpendicular magnetic field, only 
slightly changes the entropy and specific heat of the sys- 
tem. Second, we have found that the effect of barriers 
at the NS interfaces is strong and small values of barrier 
resistance, compared to the resistance of the N region, 
are necessary to observe a clear phase-dependent specific 
heat. Furthermore, assuming that the measurement is 
performed employing the method reported in Refs. j^'^i^^ 
we have taken into account the additional contributions 
to the measured specific heat in the actual experimen- 
tal setup, in which the SNS junction is embedded in a 
superconducting loop. Using Al as superconductor and 
Cu as normal-metal material, with reasonable geometry 
and sizes, we checked that a clear phase-dependent sig- 
nal is measurable. We note finally that our results for 
the electronic specific heat can be relevant for phase tun- 
ing the thermodynamic properties of mesoscopic devices 
like, for instance, superconducting single-photon sensors 
based on proximity effect.— 



IV. CONCLUSION 

In summary, in an ideal SNS mesoscopic junction both 
the electronic entropy and the specific heat exhibit a large 
deviation with respect to the normal and superconduc- 
tor case, and show a marked phase-dependence;^ In this 
paper we have studied in details such thermodynamic 
quantities in a realistic SNS junction as a function of tcm- 
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